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Abstract 

The correspondence between ordinary differential equations and Bethe ansatz 
equations for integrable lattice models in their continuum limits is gener- 
alised to vertex models related to classical simple Lie algebras. New families 
of pseudo-differential equations are proposed, and a link between specific 
generalised eigenvalue problems for these equations and the Bethe ansatz 
is deduced. The pseudo-differential operators resemble in form the Miura- 
transformed Lax operators studied in work on generalised KdV equations, 
classical W-algebras and, more recently, in the context of the geometric Lang- 
lands correspondence. Negative-dimension and boundary-condition dualities 
are also observed. 

PACS: 03.65.-Ge, ll.15.Tk, 11.25.HF, 11.55.DS. 

Keywords: conformal field theory, Bethe ansatz, pseudo-differential equations, spec- 
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1 Introduction 



A recent observation [Ij has established an unexpected Hnk between two dimen- 
sional conformal field theory (CFT) and the theory of ordinary differential equa- 
tions. This rests on a correspondence between the transfer matrix eigenvalues of 
certain integrable models (IMs), in their conformal limits [21 13], and the spectral 
determinants [U [5] of ordinary differential equations. 

The initial results [U [6l O [8] connected conformal field theories with Virasoro 
central charge c < 1 with Schrodinger problems for one-dimensional anharmonic 
oscillators. These conformal field theories are naturally associated to the Lie alge- 
bra Ai, but a generalisation to models related to An-i, with additional extended 
W-algebra symmetries, was soon established [9|ll0 t [TT | ll2]. However, contrary to 
initial expectations, a simple Lie-algebraic structure did not emerge immediately, 
and the extension of the correspondence to the theories associated with other 
simple Lie algebras g has proved surprisingly elusive. 

The purpose of this paper is to begin to fill this gap, by establishing a link 
between CFTs related to the classical simple Lie algebras and spectral problems 
associated with a set of ordinary (pseudo-) differential equations. 

We shall also prove for g = An-i, and conjecture for the other simple Lie alge- 
bras, the existence of closed systems of functional equations (V'-systems) among 
uniquely-defined solutions ip^^^ 1!)^"^^ . . . , ip^^^'^''^^^^ of a set of ranki^Q) pseudo- 
differential operators, with each pair ■i/'^") /(operator)^ being naturally associated 
to a node of the Dynkin diagram. These ■0-systems are very similar to the systems 
of functional relations introduced by Mukhin and Varchenko in the framework of 
the Bethe ansatz (B A) method for g-XXX quantum spin chains [13^ [T^ [T5t I16j , 
and in the context of the so-called Miura-opers related to the geometric Langlands 
correspondence (see, for example, [I7l[l8]). This similarity is related to the fact 
that the homogeneous 'differential' parts of the operators studied here resemble, in 
form, the Miura-transformed Lax operators introduced by Drinfel'd and Sokolov 
in their studies of generalised KdV equations and classical W-algebras [19j . 

The rest of the paper is organised as follows, ^gathers together some known, 
or easily deduced, properties of the Bethe equations for g-type quantum spin 
chains in their continuum limits. Our main results are summarised in ^ while 
extra details and numerical support for the specific A, D, B and C proposals are 
given in ^to ^respectively, and ^contains our general conclusions. There are 
two appendices: appendix |A] deals with the semiclassical analysis for ^i-related 
ODEs in the presence of string solutions, and appendix [B] describes a simple 
algorithm useful for the numerical solutions of the differential equations. The 
algorithm is a generalisation of Cheng's method from Regge pole theory |20j , and 
relies on an elegant dual formulation of the relevant boundary problems. 
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2 The Bethe Ansatz equations and their string solu- 
tions 



For any simple Lie algebra q of type An-i to (^2, a set of Bethe ansatz equations 
(BAEs), depending on a set of rank{Q) twistparameters 7={7a}, can be written 
in a universal form as [HI [22l [231 [Ml US [260 

n ^^"^^ W =-1' ^ = 0,1,2,... (2.1) 

where 

Q(")(E,7)=QW(0'=i?,7) , (2.2) 

and the numbers -E^'*^ are the - in general complex - zeros of the functions Q^"'^: 

QW(i^W,7)=0. (2.3) 

In (j2.ip and (j2.3p the indices a and h label the simple roots of the Lie algebra, 
the matrix Bab is defined by 

Bai,= , ^""'"'j 12 , a,6 = l,2,...,ranM0) (2.4) 
I long roots I ^ 

and the a's are the simple roots of g. The constant J7 is a pure phase, parame- 
terised in terms of a real number ^>0 as 

= „p(.|L) (2.5) 

where hy is the dual Coxeter number of g. Strictly speaking, the BAE ()2.ip arise 
from taking a suitable continuum or field theory limit of the lattice model BAE, 
in the fashion explained in, for example, [27]. 

The functions Q^"^ appearing in (j2.ip have a characteristic asymptotic be- 
haviour at large values ol —E 

lnQ(")(-j^,7)=m. .'^f^^ A -EY + ■ ■ ■ ■ (2.6) 

For the Cn and models the sets {ma} are given in table[lll The 

only free parameter -the overall constant m in table [1]- depends on the way the 
conformal field theory limit is reached. 



* For finite lattice models, the explicit diagonalisation of the cases has been performed 

through the algebraic Bethe ansatz by Schulz [21] and also by Babelon, de Vega and Viallet |22] . 
For Cn and Dn models, it has been done by Reshetikhin [231 124) . There is a shortcut to reach 
the same conclusions via the so-called analytic Bethe ansatz of Reshetikhin [25) . and Wiegmann 
and Reshetikhin [26J. 

^ The constants {ma} are related to a particular matrix Kab emerging from the analysis 
of the Bethe ansatz. For simply-laced algebras, Kab is proportional to the Cartan matrix and 
v={mi,m2, . . . , rrir) is its Perron-Frobenius eigenvector. 
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Model 




ma 


An-l 


n 


rua = 2m sm{a-^) , (a = 1, . . . , n — 1) 


Dn 


2n-2 


m„_i = m„ = m, nia = 2msm{a-^), (a = 1, . . . , n — 2) 


Bn 


2n - 1 


nin = m, rua = 2m sm{a-^), (a = 1, . . . , n — 1) 


Cn 


n + 1 


ma = 2msm{a2^), (a = l,...,n) 



Table 1: Dual Coxeter numbers and coefficients {ma} for models based on classical simple Lie 
algebras. 

The negative real E axis is also the direction of maximal growth for In Q^""^ (E) 
as \E\ — > oo. Prom (j2.6p . the Hadamard order of Q^"^ is therefore /i and, in the 
so-called 'semiclassical regime' < fi < 1, Q^""^ can be written as a convergent 
infinite product over its zeros as 

Qi-){E,^) = Q(^H0,j)f[(^l-^^ . (2.7) 

It turns out that the Bethe ansatz roots generally split into multiplets with ap- 
proximately equal modulus and that the ground state of the quantum spin 
chain corresponds to a 'pure' configuration of roots, containing only multiplets 
with a common dimension 

4 = 7^. (2.8) 

Dan 



The integer K in (12. Sp depends on the particular spin chain under discussion, 
and corresponds to the degree of fusion [281 EHl E^- Foi' 5=^i ) the spin-j A\ 
quantum chains, K=d\=2j . 

It is generally expected [31] that for large values of i the zeros asymptotically 
tend to the perfect string configurations: 

argi?f)~(4 + l-20|^ , I = 1,2,..., da. (2.9) 

Appendix A contains some further discussion of the asymptotic behaviour of these 
string solutions. 
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3 Summary of the main results 



This paper is about the correspondence between the Bethe ansatz equations (12. ip 
for g=An-i, Bn,Cn, Dn and spectral problems associated to solutions ip^XjEjg) 
of particular pseudo-differential equations with vanishing boundary conditions 



g) = o e 



(M > K/{h^-K)) 



(3.10) 



imposed at large x on the positive real axis. To specify these equations we intro- 
duce the n*'^-order differential operator [11] 



n-l 



{n-l))D{gn-2 

D{g) = 



-(n-2)) ... D{gr-l)D{go) 

_d_ _ 5 

dx X 



(3.11) 
(3.12) 



with 

g={5n-i, • • • ,5i,5o} , g"^={?^ - 1 - 50,?^ - 1 - 51,' 



,n 



1 - gn-i} . (3.13) 



We also use the inverse differential operator (d/dx) ^, generally defined through 
its formal action 



d 

dx 



„s+l 



S + l 



The following properties hold 



d 

dx 



d 

dx 



X 



d 

dx 



dx 



X 



(s + 0), 



and the integration by parts property 

-1 



d 
dx 



dx 



f{x)g{x) - ( — 



g{x)—f{x) 

dx 



(3.14) 



(3.15) 



(3.16) 



is satisfied, provided in the x-expansion of f{x)g{x) about the origin the constant 
term is absent. In the following we shall assume the validity of (j3.16p by working 
implicitly with non- integer values for the parameters gi introduced in (j3.13p . and 
by invoking continuity of the final results in these parameters. 
Finally, we define a basic 'potential' 



Pk{E,x) = {x 



h'^M/K 



E) 



K 



(3.17) 
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With this notation in place, the following pseudo-differential equations are 
the main concern of this article: 



An-i (su(n)): 

■(-irZ)„(g) - PK{x,E)ji;{x,E,g) = , (3.18) 
with the constraint Yl^=o 9i= ^^\^^ ■ 

Dn (so(2n)): 

(^Dn{g^) (^^) Z?„(g) - ^Pk{x,E) {J^ ^Pk{x,E)^ i;{x, E,g)=0 

(3.19) 

Bn (so(2n+l)): 

(Dnig^)Dnig) + VPk{x,E) (^-^^ VPk{x,E)^ i;{x, E,g)=0 (3.20) 
Cn (sp(2n)): 

(^n(g^) Pk{x,E)\ i;{x,E,g)=0 (3.21) 



The correspondence we propose links the ground-state Q^^^'s of (j2.ip to par- 
ticular solutions ip (|3.10p of equations (|3.18H3.2T]) . 

In order to clarify this statement we introduce an alternative basis of solutions 
{Xi{x, E,g)} to (|3.181l3r2T]) . characterised by their behaviour near the origin 

Xiix,E,g) ^x^^ + ... , x^O, (3.22) 

where the A's are the ordered (Aq < Ai < . . . ) roots of the appropriate indicial 
equation (see table [21). 

Writing -0 as a linear combination of the x's, we have in general 

^{x,E,g) =Q'^l{E,g)xoix,E,g) +Q\ll{E,g)xi{x,E,g) + ... . (3.23) 

If the zeros of Q^^^ {E, g) are {E^^^}, then for Eg{eI^^} the function x^^'^ip{x, E, g) 
vanishes exceptionally both at x=oo and at x=0. This allows the coefficient 
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Model 


indicial equation 


^n— 1 


nr="o(A-5.) = o 


Dn 




Bn 


nr=o(A-5.)(A-/.^+5,.) = o 


Cn 





Table 2: Indicial equations. 



Qjgj^(£', g) to be identified with the spectral determinant for a boundary problem 
defined on the positive real axis (see, for example, [ll[5]). An alternative (dual) 
definition of the spectral functions Q^^{E,g) in terms of the adjoint equations 
to (|3.18|) - (j3.21|) is briefly discussed in appendix [Bl 

We claim that for classical Lie algebras with arbitrary degree of fusion K, the 
ground-state Q^^\E,^ys in (|2.ip and the functions Q^^^{E,g) in ()3.23p coincide 
up to a trivial normalisation, so that 



(3.24) 



Q^^HO'7) Q«(0,g) 
Moreover, from the WKB approximation 

lnQ\^^{-E,g)=K{-Ef + ... (^>0) (3.25) 

with 



_ K(M + l) fh^M , , 

u = — -, k = k[ , — (3.26) 



and 



^ ' Jq V ' J r(l + l/o + l/6)sin(7r/6 + 7r/a) ^ ' 

Therefore, in order to have a match between (j3.25p and (j2.6p . we must set 

sin(^5ii) ^ __ ( ^ 2ttM 



Given a particular ordering convention, the relationship between the twist pa- 
rameters {7a} and the constants {ga} is given in table [3l 
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Model 


ordering \/ i < j 


{9a} ^ {la} 


A 


9i < 9j 


7a — " yl^i=0 9i 2 J 






7a — a \ 2_ji=Q 9i — 2'^ j ) l^a — i, . . . , n — 


Dn 


Qi < gj < /2 


7n-l = 2 [l^i=0 9i-2^ ) 

^ = " f y"-2 g. _ g„ 1 _ n^h"^] 
2 \ Z--/2=U t/n—i 2 } 


Bfi 


Qi < Qj < /l^/2 


7a = «(E?=o5.-f/i^) 


Cn 


gi < Qj <n 


7a = a (E^=o 9i - anj 
7n = f (Er=-o<7.-n2) 



Table 3: The relationship between the set of parameters {ga} ^ {la} with a=2K/Mhy . 

Various consistency checks, including the WKB approach and numerical work, 
support the correspondence both qualitatively and quantitatively. 

Finally, starting from equations (|3.18p to (j3.2ip . the Bethe ansatz equa- 
tions and table [3] were obtained with the help of a system of functional rela- 
tions involving -E, g) = 'ip{x,E,g) together with other auxiliary functions 
ip'-''\x,E,g) , (a = 2,...,rank{0)) (see Sl-g7]and PE]). We sell 

^pi^^ =i;'-''\u;^x,n''E,g) (3.29) 

where 

/ 2irM \ 

as in (j3.28p . and 

For the simply-laced algebras the ^—systems can then be written in the compact 
form 

rank(Q) 

W[i^H,^['''']= n (^^'^)^" ' (3.32) 

2 2 



•'■In the An-1 models ipl^\x,E,g) = a;'" ^^'''/^i/_fc(a;, _E, g), where is the function defined 
in §3 of [11]. 
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where Aab='^Sab — "^Bab is the incidence matrix of the corresponding Dynkin dia- 
gram and W the Wronskian: 

W[f,9] = f{x)^g{x)-g{x)^f{x) . (3.33) 
ax ax 

Equation (j3.32p is proven in ^4.2l for 0=A„_i, and our numerical resuhs indirectly 
support the validity of (|3.32p for g=Dn- Currently, we have no analogous pseudo- 
differential equations for the exceptional Lie algebras but the similarity between 
(j3.32p . the relations proposed in [I3l [Ml [TSj and the other functional equations 
(Y-systems and T-systems) discovered in the framework of integrable models 
[32l [33l [34l [35l [36l [37] suggests the validity of (|3.32p in its more general form. For 
the non simply-laced algebras our conjectures are 



Bn : =^('^-i)v('^+i); a = l,...,n-l, 

2 2 

w^ivil , ^^r^] = viv'V^r^'^ • (3.34) 

4 4 4 4 

(Our root convention is (aj|aj)=2 for i = 1,2, ... ,n — 1 and (a„|a„)=l.) 



4 4 

4 4 4 4 

2 2 

(Here, (aj|ai)=l for i = 1, 2, . . . , n — 1 and (an|an)=2.) 



.,n-2. 



(3.35) 



4 4 4 4 

2 2 
2 2 

(Here, (ai, ai)=(a2> 02)=! and (03, a3)=(a4, a4)=2.) 



(3.36) 



G2 : Ty[vi'),vl'^]=V(2) ^ 

2 2 

6 6 6 6 

(Here, (ai|ai)=3 and (a2|a2)=l-) 



(3.37) 
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Again, these relations are not proven but we have indirect numerical evidence 
for Q=Bn, Cn- Further details and numerical support for the above conjectures 
are provided in the following sections, which examine the A, D, B and C cases 
in turn. 



4 The An^i models 

The ODE for the An-i models is 

-l)"+iD„(g) + Pk{x, E) ) i?, g) = , (4.1) 



where the operator Dn{g) and the generalised potential Pk{x, E) were defined in 
^ and the additional constraint 



^ nin—l 



i=0 

ensures that the term dx'^-i is absent. 

The function 'ip{x,E,g) is defined to be the most subdominant solution on 
the positive real axis, with asymptotic behaviour, for M > K/{h'^ —K), given by 

i>{x, g) ~ AA x(^-")^^/2 exp{-x^^+'^/{M+l)) (4.3) 

as a; ^ cxd. The coefficient M represents an E- and g-independent normalisation 
constant. 

The K=l cases have been extensively discussed in [9l [lOl [TTl [12]; they are 
related to the WAn-i conformal field theories with integer Virasoro central charge 
c = n— 1. Alternatively, at particular values of the parameters g and M they can 
also be put in correspondence with the minimal coset conformal field theories 

(In-l)l X (In-i)L ^^^^ 
iAn-l)L+l 

with a simple relationship between L in (j4.4p and the parameter M in (j3.17p . 

The generalisation to integer K>1 comes from an observation by Sergei Lukyanov 
|38j for the Ai case, for which numerical and analytic support was later provided 
in [39] and in [30]. It is reasonable to conjecture that this generalisation works 
both for An-i with n>2 and, up to minor modifications, for the other models to 
be discussed in this paper. In analogy to (j4.4p . at particular values of g and M, 
the integer K>1 cases should correspond to the cosets 

, (4.5) 

Qk+l 

which describe conformal field theories with Q=An-i, B^, Cn, D^- In appendix [Al 
we shall explain, in the simplest case, why potentials of such forms naturally lead 
to string patterns of roots of the sort mentioned in the introduction. 
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4.1 Negative-dimension dualities 



It is interesting to note that there are formal duality relations among our pseudo- 
differential equations involving negative values of n and K. Consider the An-i 
ODEs with the twists g={0, 1, . . . , n — 1}: 

((-1)"+'^ + PK{x,E))i^(,x,E) = . (4.6) 

Setting '4){x,E) = Pk{x, E)il;{x, E) and multiplying from the left by (^)~'^, the 
result is 



(7 \ — n 
+P-Kix,E))ti^ix,E)=0 



(4.7) 



Comparing (j4.7p with (j4.6p and taking into account the boundary conditions, we 
see that there is a formal duality and a spectral equivalence between the initial 
n*^-order ODE and the pseudo-differential equation (j4.7p : 

{n,M,K} ^ {-n,M,-K} . (4.8) 

Though the above manipulation might look purely formal, it strongly resembles 
previously-observed W-algebra dualities |i41j: 

(l_„)x X (1_„)l iAn)-K X (A„)_i 



(^-n)i<'+L {A-n)-K-L 

A discussion of the precise relation between L and the ODE parameters {n, M, K, g} 
is not important for the current naive considerations and we shall postpone it to 
the future. 

Whilst the duality ()4.8p remains at the moment a purely formal observation, 
the second duality discussed in [H] 



{D-n)K X {D-n)L _ {Cn)-K/2 X { Cn)-L/2 



{D~n)K+L (Cn)- 



K/2-L/2 



will lead us to the proposal (j3.2ip for the C^-related equations. (For the quan- 
tisation of the classical Dn W-algebras in relation to the Miura-opers see, for 
example, ^2].) 

Negative-dimension dualities resembling those described here are also well- 
known to group theorists; many more details can be found in 



4.2 Auxiliary functions and the ^/'-system 

The solution ^(x, E, g) of ()4.ip with the asymptotic behavior (j4.3p is not the only 
function associated to the A^-i Bethe ansatz equations. In order to derive the 
Bethe ansatz itself a total of n—1 functions • • • > one for each node of 
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the An-i Dynkin diagram, should be introduced. These auxiliary functions are 
solutions of generally more-complicated ordinary differential equations. 

Following [11], all the functions {x , E , g) decaying at large x on the posi- 
tive real axis can be constructed from ^ = ijj^^^ as 

^(a) ^ p^(a)[1^^3_a^^^^^a_l] ^ [-0 , "03^ , . . . , t/^o^] , (4.11) 

where a = 1,2, ... , n—1, W^""^ [fi, . . . , fa] denotes the generalised Wronskian of 
the set of functions {fa} 



W^''\fu...,f a] =det 



a-l 



f,-f,...,^f 



(4.12) 



with / = (/i,/2,...,/a) (so that W^^\f,g\ = W[f,g], cf. eq.dSjH])), and 
denotes the 'rotated' solution (|3.29p . 
Finally, normalising by 



iin-l)/2 



M= ^ =^ iP^''\x,E,g) = 1 . (4.13) 



n 



Since is a solution of an n*'^-order ODE, a naive counting of degrees of 
freedom shows that the order of the ODE satisfied by ip^""^ should be 

n' 

(4.14) 



(n — a)l a! 

The (n—a+l) and the (a) equations are related by a (g ^ gt)-conjugation [51 [TT]. 
arising from the Z2 symmetry of the Dynkin diagram. Notice that the result (j4.14p 
exactly matches the dimensions of the basic representations of An-i; these are 
again in one-to-one correspondence with the nodes of the Dynkin diagram. 

Fortunately, in order to derive the Bethe ansatz equations an explicit knowl- 
edge of the remaining n—2 ODEs is unnecessary: the derivation of [llj was instead 
based on the Stokes relation associated to (14.11) 



^i-l)^C^>'\E,g)ykix,E,g)=0, (4.15) 

k=0 

where, according to [TT] , 

yk{x, E, g) = a;("-i)^/2v.(a;-'=x, Q-'^E, g) , (4.16) 

C^^\E,g) = 1 and the Stokes multipliers C^'^\E,g) with /oO are analytic func- 
tions of E and g. Stokes relations for the D„, -B„ and C„ equations (|3.19p . (j3.20p 
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and (j3.2ip also exist, but we have encountered some subtle complication^^ in 
generalising the approach of [H] to these cases. 

Instead, the strategy here is based on the conjectured validity of a simple 
'universal' system of functional equations among the ip^""^ functions, which leads 
immediately to the Bethe ansatz equations and bypasses the analysis of Stokes 
relations. 

We shall now prove the An-i ■i/'-system (|3.32p : 

^(a~l)^(a+l) ^ ^[^(a) ^ ^(a)j ^ ^(Q) ^ ^(n) ^ ^ ^ (^^-^^^ 
2 2 

The proof is based on the observation that determinants satisfy functional equa- 
tions, in particular the so-called Jacobi identity 

A A[p, q\p, q] = A[p\p] A[q\q] - A[p\q] A[q\p] . (4.18) 

Here, A is the determinant of an {a + 1) x (a + 1) matrix and A\pi,p2\qi, q2] 
denotes the determinant of the same matrix with the pi^2 — th rows and qi^2 — th 
columns removed. 

In order to prove ()4.17p we have to consider three different cases: a=l, 
l<a<n—l and a=n—l. The a=l case follows from the definition of ■0^^^ given 
in (14. lip . Equation ()4.17p for l<a<n—l follows from the following chain of iden- 
tities 

b 

_( 1 'i(a-l) TT/(a+l) r a-2 a- 2 _a ai TT^(a-l) r a-2 a— 2 i 

— y i.) vv [ 2)---)2' 2'' 2\ L 2'''''2J 

= (_l)(«-i)AA[a,a + l|a,a + 1] (4.19) 
where we have identified 

A = VK(«+i)[-^,...,«^,-|,|] (4.20) 

and 

A[a,a + l|a,a + l] = [-«^, . . . , (4.21) 

This is nothing but the LHS of the Jacobi identity (|4.18p . Then an application 

^In the D„ case these comphcations are probably a consequence of the fact that the ODEs 
associated to the Z2-conjugate nodes in the Dynkin diagram are somehow more fundamental 
than eq. (|3.19|l . This latter equation is more naturally associated to the first node on the 'tail' 
of the diagram. 
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of the Jacobi identity naturally proves ()4.17p in the following way: 

-Q(^W)A,, ^ (_i)(«-i)(A[a|a] A[a + l\a + 1] - A[a\a + 1] A[a + l\a]) 



— r '^"2 g— 2 a] t;[/(i) [ a— 2 a— 2 an"! 

[ 2 5 • • • ) 2 ' 2J '^''^ ^ 2 ' ■ ■ ■ ' 2 ' 2V 

= . . . , f ] VF('')[-f , . . . , ^] 

TT/'{i) r ^1 a— 2 g— 2 1 T^/{g) [ g— 2 a— 2 a] 

I'l' L 2 ' 2 ' • • • ' 2 J ''^'^ ^ 2 ' • • • ' 2 ' 2J 

= _ ^'(-)^[-) = ^[^W , ^W] . (4.22) 

22 22 22 

Finally, for a=n—l the previous calculation shows that 

^(n-2)^(n) ^ p^j^(n-l)^^C-l)] , (4,23) 
2 2 

Choosing M = '^^"^^^ gives ■0(")=1 and (j4.17p is proved. 
4.3 The Bethe ansatz equations 

In this section we shall show that the BAEs are a simple consequence of the 
■0-system. First, recall the alternative x-basis of solutions ()3.22p and the formal 
ordering of table [3] 

gi<gj y i<j . (4.24) 
These solutions are defined by their behaviour as 2; ^ 

X.(x,^,g)~x^'+0(x^'+") , Xi = gi . (4.25) 
Next, expand ip{x,E,g) in the x-basis 

n-l 

ij{x,E,g) =Y,Q\l]{E,s)x^{x,E,g) , (4.26) 

i=0 

and use the property 

,k„ r\k TT' „\ , ,k\ 



Xiiuj'x, n'E, g) = oj'^'-^Xiix, E, g) , (4.27) 
to obtain 

n-l 

Mx,E,g) = J2Q''ll]i^'E,g)u;''>^'x^ix,E,g) . (4.28) 
Now expanding the determinants in the determinants for this new basis leads to 



''>{x,E,g) 



T.'Q[:L g)^^"^ fc^ , . . . , X.^ ] (4.29) 
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where j = ^-g^, . . . , Yli denotes the sum from to n— 1 of the set {ij} 
while in there is the additional constraint < ii^ < < • • • < ^a-i • 

_ 2 2 2 

A family of x-independent equations is obtained by identifying from the LHS 
and RHS of (j4.17p the terms corresponding to the same power. It is possible a- 
priori to identify from every determinant only the highest, second highest, lowest 
and second lowest orders. We shall extract the leading orders, though similar 
results can be obtained from the subdominant ones. 

Setting 

Qt\E,g) = Q\^l^^^^^{n'^E,g) , Q^:\E,g) = g[oi.,„_2,a](f^'^>g) (4-30) 
we have 

i;i^\x,E,g) = ^^"«Ql")(i?,g)Ty('^)[xo,...,Xa-i] + 

^fcK+i-A._i)gW (^^ g)^(a) ^^^^ ^ Xa] + ... . 

(4.31) 

The orders of the first and the second terms in (j4.3ip are a^— a(a— 1)/2 and 
aa+i—Xa-i—a{a—l)/2 respectively, where 



a-l 
i=0 



a = Y,Xi. (4.32) 



Substituting ()4.3ip into ()4.17p . comparing the leading terms of both sides for 
small X and using the relation 

(also proved through the Jacobi identity), we find 

g(a+l)g(a-l) ^ ^i(A.-A._i)g(a) g(a) _ ^|(A._i-A.)g(a)g(a) ^^^33^ 

22 22 

g(0) = 1. Finally, let eI""^ be a zero of Q^""^ {E, g). Evaluating the above equation 
at E = ^l^/'^E^"'^ and at E = Q.^^/'^E^°'\ and by dividing the two equations thus 
obtained, we find the An-i Bethe Ansatz equations 

n^^-^ ^w =-1- .=0,1,2,..., (4.34) 

ab 



b=l Q-B.S^: 



where 

2K ("^ (n-1) 

la- " " 



a=0 
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Notice that in writing relation (j4.35p we have used the identity 



^{X,-X,_,) = -Y^Bablb, (4.36) 

b=l 

the ordering (j4.24p and imposed the constraint 7o=7n=0. 

As shown in table [H for K=l there is very good agreement between the 
IM results obtained from the solution of a suitable non-linear integral equa- 
tion (NLIE) (see §6 in [llj) and the direct numerical solution of the ODE. 



Level 


A4 NLIE 


ODE numerics 


^0 


14.0495626907 


14.0495626922 




47.7146839354 


47.7146839363 




95.1785845453 


95.1785845456 




154.202021469 


154.202021470 




223.483044292 


223.483044292 



Table 4: Comparison of IM results {A4 NLIE) with the direct numerical solution of the A4 
ODE with K^l, M=10/21 and g=(0.2, 1.02,2.3,3.421). 



Table [5] also shows the good agreement between the IM results obtained using 
the spin-1 NLIEs OSS] and the exact solution of the ODE with K=2 (see ^Ol ). 



Level 


Ai NLIE with K=2 


ODE numerics 


^0 

E^^ 

E^^ 


1.49259741085 ± 1.60304589242i 
2.31180377628 ± 2.38537059826i 
2.91183770898 ± 2.97068128676i 
3.40837129214 ± 3.45880577384i 
3.84143464742 ± 3.88626414305i 


1.49259741085 ± 1.60304589242i 
2.31180377628 ± 2.38537059826i 
2.91183770898 ± 2.97068128676i 
3.40837129216 ± 3.45880577388i 
3.84143464640 ± 3.88626414641i 



Table 5: Comparison of IM results {A\ NLIE) with the exact solution of the A\ ODE with 
K=2, M=l and go = 0. The set {2S*^'} are the exact eigenvalues of the Bi Unear ODE of gOI 



5 The Dn models 

The Dn pseudo-differential equation ()3.19p is 



^n(gt) 'l?n(g) -^/Pk{x,E) ^Pk{x,E)^ ^{x,E,^) = . 



(5.1) 
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Following the An-i example, we start from the solution ip{x,E,g) of (jS.ip with 
asymptotic behaviour 

V'(x,£;,g)~AAx-^ *^/2exp(^-^— , {M>K/{K'-K)) (5.2) 

as X — > oo on the positive real axis, and introduce the alternative basis of solu- 
tions ([3:22]) 

X.(x,^,g) , i = 0,l,...,2n-l (5.3) 



characterised by their behaviour near the origin 

X.(x,i?,g)~X^'+0(x^'+2n-l) ^ 



Xi= gi , for i < n - 1, 

\i = h>' - g2n-i-i , for i > n - 1. 

(5.4) 

In ()5.4p the parameters Aj represent the 2n solutions of the indicial equation 
(see table [2]) with the ordering 

gi < 9j < /i^/2 , \i<Xj , y i<j . (5.5) 

5.1 The ^/'-system and the D„ Bethe ansatz equations 

To extract the D„ BAE, we start from 

V-i'^ = ^fe(x, E, g) = ^{iv^x, n'^E, g) , (5.6) 

where Q and uj are as defined in (j3.28p and (|3.29p . and assume the validity, for a 
suitable value of the normalisation constant M, of the ^-system (|3.32p 

= a=l,...,n-3 

2 2 

T^lV'iV'^V'l""'^] = V'^^-'^V'^""'^^"^ (5.7) 
2 2 

22 22 

Equations (|5.8p and Jacobi identity (jS.Sp used in reverse imply the following 
relations linking the remaining functions 'ip^"'\x , E , g) to ip^^\x, E,g): 

^(a) ^^(a) ^^{a)r^^^^^^___^^^i a=l,...,n-2, (5.8) 

2 2 2 

^(r^-i) =^(«-i)^(") = |y("-^)[V'2^,V'4^,...,V'^] , (5.9) 

2 2 2 

and 

^(n) ^^(n-l)^(n-l)^^H^{n) ^^(n)[^^^^^^^^^^^^] . (5.10) 
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Now notice that the auxihary functions (j)^"'\x , E , ^) defined in (|5.8p . (j5.9p and 
(jS.lOp satisfy an ^-type V'-system 

^(a-l)^(a+l) ^ ^j^(a) ^ ^(a)^ ^ a = 1, . . . , n - 1. (5.11) 

2 2 

Therefore, the arguments appHed in ^4.31 go through in the same way: 



{E, g)W^-'> [xo, . . . , Xa-2, Xa] + ... (5.12) 



(a = 1,2,... ,n — 1). The orders of the first and the second terms in (j5.12p 
are given by Oq— a(a— 1)/2 and aa+i— Aa_i— a(a— 1)/2 respectively, with = 
E"=o and 

Qi")(i?,g) = g('^)(17'^i?,g) , Q^^\E,g) = Q(-\n''E,g) (5.13) 

are (/)-related spectral determinants. 

Using equation (jS.lip we establish the following identity among (/)-related 
spectral determinants 

q('*+i)(e)q('^-i)(e) =a;^(^"-^«-i)Q5!i(^)Qi"^(^)-w^^^""'"^"^Qi"^(^)0il(^); 

2 2 2 2 

(5.14) 

which leads to 

z.]: ]: = -ot(^'"-^'"-) (5.15) 

g(a-l)(^(a)^ gW^^Ja)^ ^(^^^(i^f )) 

2 2 

with a = l,2,...,n — 1 and a= . We then make the following identifications 

Q(o)(i?,g) = Q^''\E,g) = \- (5.16) 

Q('^)(i?,g) = Q(")(^,g), (a = l,...,n-l); (5.17) 

Q("-i)(i?,g) = Q("-i)(i?,g)Q(")(i?,g) ; (5.18) 

g(")(i?,g) = Q(T^)(i?,g)Q(r'^(i?,g), (5.19) 

2 2 

which reflect the relations among the (/)'s and the -(/''s in the above. In (|5.19p we 
have implicitly assumed 

V'(")(x,ii;,g) =o(V'("-^)(x,i?,g)) (5.20) 

— > and also that (see the discussion in ^5.2p 

Q("-l)(£;, {ffo, . . . , 5n-2, <?n-l}) = Q^"^^, {50, • • • , 9n-2, /l'' " 5n-l}) ■ (5.21) 
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Plugging relations (|5.17p . (|5.18p and (|5.19p into (|5.15p and using (|5.2ip is it easy 
to check that (I5.15P can be recast in the universal form ()2.ip 



n^^"" ^(6) A =-^' . = 0,1,2,... (5.22) 



where is the D^-related matrix defined according to (12. 4p and we have im- 
posed the extra condition 

7n - 7n-i = n - 1 - (5.23) 

to fix the exact {ga} <-> {7a} relation as given in table [3l 

Condition (I5.23|) guarantees that when gn-i—n+l=0 the operator (d/dx)^^ 
in (|5.ip acts directly on a d/dx and the relevant equation reduces to an (2n— 1)- 
order ODE. When this occurs 7n=7n-i and so the pair of Z2-conjugate nodes of 
the Dynkin diagram are 'twisted' in exactly the same way. Further, a change of 
sign in the RHS of (j5.23p swaps 7„ and 7n-i, a property that naturally reflects the 
presence of the Z2-symmetry in the Dn Dynkin diagram. All these properties are 
confirmed by the analysis of §5.2l and §5.3l and by 12-digits of numerical agreement 
at K=l with g={0, 1, . . . , n - 1}, 

/(a-l)a /i^N n/i^ , 
7a = a I ^ a— 1 , 7n = 7n-i = ^ , (5.24) 

between NLIE and ODE results. 

Table[6]shows the (still-excellent) agreement at K=l away from the 7„ = 7n-i 
surface. Appropriate K>1 NLIEs are unknown but numerical results qualita- 
tively reproduce the expected IM scenario of ^ Further analytic support to the 
proposed ODE/IM correspondence for Dn is given in ^5.2^ ^5.31 and ^5.41 below. 



Level 


D4 NLIE 


ODE numerics 


^0 


17.8625636061 


17.8625636061 




50.2942213433 


50.2942213430 




92.8267466445 


92.8267466442 




143.348705065 


143.348705065 




200.738324171 


200.738324172 



Table 6: Comparison of IM results {D4 NLIE) with the direct numerical solution of the D4 
pseudo-differential equation with K=l, Ai—1/3 and g — (0.2, 1.1, 2.3, 2.95). 
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5.2 Example 1: 1^2 ~ © 

The D2 algebra can be decomposed into a pair of independent Ai algebras, mir- 
roring an analogous factorisation in the BAE. In this section we shall prove that 
the solution il){x, E, g) to (jS.ip with n=2 is the product of two solutions of Ai- 
related ODEs. 

We start from the general D2 equation: 

A. + ^\ f ^ . 9i-'^\ f f d gi-l\ { d go 



dx X J \dx X J \dx J \dx x J \dx x 

-^/Pk{x, E) VPk{x, E)^ ij{x, E,g)=0. (5.25) 

Expanding and integrating by parts, we obtain an equivalent equation: 
d^ d ^dP 



f d: ^ d dP , ^ 



where for this subsection it is convenient to define P{x, E, k) = ^[{x^^'^ / ^ — E)^ + 
^] , and g' =gl- 2go + gf - 2gi + 1. 

We now set Z{x) = Xi{x)X2{x), a product of the solutions of two Ai (spin-^) 
equations, which for general p and a satisfy 

-^Xi{x,E,p) = P{x,E,p)xiix,E,p) , (5.27) 

—X2{x,E,a) = P{x,E,a)x2{x,E,a) . (5.28) 

To show that Z{x) satisfies equation (|5.26|) . we differentiate and repeatedly use 
the Ai equations to find 

d^Z dP , „ p + cr^ „ p + a, dZ „ ^ dx2 

+2P{x,E,a)^X2 ■ (5.29) 
ax 



If we now define the Wronskian 



dX2 dx 

vv = XI 

and use 



.rr ax2 axi on^ 
W = xi-^ -1-X2 (5.30) 

ax dx 



dX2 I ( dZ \ dxi I (dZ \ 
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then ()5.29p can be written as 

d^Z dP p + a p + a dZ p-a 

dx'^ dx 2 2 dx 

In order to express W in terms of Z we differentiate, apply ()5.27p and (I5.28p . and 
then integrate: 

dW _^^d\2 d^Xi^^_cT-p^ ^ w = (a-p)(—y'— (5 33) 
dx dx"^ dx"^ \dx J 4x^ 

The resulting equation for Z, 

d^ . AT^/ 1^ P + '^\d ^dP , ^ p + a. 



16a;^ \dx J 

exactly matches equation ()5.26p provided the following relations between qq and 
gi, and /? and a hold: 

p + a = 2{gl - 2go + gj - 2gi + 1) (5.35) 
^ = (50-l)(5i-l) . (5.36) 

If p=a then either go or (71 has to be zero and the integral operator in the D2 
equation acts on a total derivative. This observation agrees with the discussion in 
^5.11 about the relation between {d/dx)~^ and an asymmetric choice of the twists 
7n and 7„_i. 

5.3 Example 2: D3 ~ 

The BAE for and D3 are the same under identification of Bethe roots. It is 
therefore interesting to discuss the exact correspondence between the two models 
at the level of the pseudo-differential equations. Actually, it was the study of this 
case that lead us to the general Z)„-related equations. 

We start from the observation that the solution of the ^3-related ODE 

D^{g)-PK{x,E)^^P{x,E,g)=0 , (5.37) 

is associated to the first node of the ^43 Dynkin diagram, while the solution of 
the D2, equation 

(^D;{g^) (J-^ ' D;{g) - WPk{x,E) (J-^ VPk{x,E)^ cl^ix, E,g)=0 

(5.38) 
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is more naturally associated to the central node of the L>3=yl3 Dynkin diagram. 
(A constant factor r has been included in eq. ()5.38p to take into account the 
possibly-different normalisations for the E parameters.) 

Therefore we are looking for a relationship between (j){x,E,g) and 

V^(2)(^,^,g) =W^[^(i),^(i)]. (5.39) 

2 2 

For simplicity, we perform the calculation at g={0, 1,2}. We set 

^(2)(a;^^) = [0,1] (5.40) 
where we have introduced the shorthand notation 

so that 

^[i,j] = [i + l,j] + [t,j + l] , [i,i]=0, (5.42) 

and 

—^|J^^^,{x,E) = -Pk{x,E)^1\{x,E) . (5.43) 
Taking derivatives five times and using the above equation we have 

= [0,1], ^V(^) = [0,2], ^V(^) = [l,2] + [0,3], 
^V'^^) ^ 2[1,3] + [0,4] = 2[1,3]-Px[0,0] = 2[l,3] , 
^V^') = 2[2, 3] + 2[1, 4] = 2[2, 3] + 2Px^(2) ^ 

|,^(^) = 2[2,4]+2A(p,^(2)) = 2P.A^(^)+2A(p,^(2)),(5.44) 
and therefore obtain the desired ODE 

- 2^/Pk{x,E)^ ^Pk{x,E)^ ^(2) E)=0, (5.45) 

which matches ()5.38p at r=2. 

We have also checked that the solution ^jJ^^\x,E,g) of the more general As- 
related differential equation (j5.37p leads to a function V'^^^ {x, E, g) = , V'l^^] , 

~ 2 2 

which is the solution of 

(osig^) ' D-M - 2^Pk{x,E) (J-^ VPk{x,E)^ V^'^ {x, E,g)=0. 

(5.46) 
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As already seen in H5.2\ in order to recast the resulting equation in the fac- 
torised form (15.460 one has perform a number of integrations by parts. 
The exact relation between the ^3 and sets of parameters is 

2go = 1 + go + 91- 92 , 

291 = 1 + 50 - 5"i + 92 , (5.47) 
2^2 = 1 - 50 + 91 + 92 ■ 

5.4 Relationship with the sine-Gordon model 

The reader may have noticed that the sets of numbers {nia} for the Dn and 
models summarised in table [1] match the mass spectra of the sine-Gordon model 
at particular values of the coupling constant. From the sine-Gordon point of 
view the Dri-related mass spectrum emerges at the reflectionless points where the 
scattering between the solitons becomes purely diagonal. This link between the 
sine-Gordon model, affine Toda field theories and perturbed coset CFTs has been 
discussed in various places jlHl HTl HU 09] . 

In this section we would like to point out that there is a simple connection 
between equation ()5.1|) taken at 

K = l, M = l/(n-l), g = {0,l,...,n-l} (5.48) 

and the Schrodinger problem associated, through the first instance of the ODE/IM 
correspondence [T, '6], to the CFT limit of the sine-Gordon model. We start from 
()5.ip with parameters ()5.48p : 

and require xi^^ to be absolutely integrable on the full real line; this restricts 
the possible values taken by to a discrete set. Fourier transforming (j5.49p yields 

("^ " 11 + ((-1)"^'"-' - ^)) X{k, E) = 0, (5.50) 

and replacing 

k^ik , E^-E , x{k, E) k-^/^x{k, E) (5.51) 

we finally find 

Equation (j5.52p exactly matches the ODE associated in [U |6] to the refiectionless 
points of the untwisted sine-Gordon model at its c=Cefr=l conformal point. 

This simple observation gives extra support to the correctness of the pro- 
posal (|5.ip . and it leads naturally to the Bn proposals discussed in the next 
section. 
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6 The Bn models 

The discussion in §5.41 of the Unk between the sine-Gordon and Dn scattering 
theories at specific values of the parameters can be extended to the Bn models [HI 
149] ■ This and further considerations led us to the ODE ()3.20p . which we repeat 
here: 

(^Dnig^)Dnig) + VPKix,E) (^-^^ ^/Pk{x,E)^ ij{x, E,g)=0. (6.1) 

The results of ^5.11 and ^5.21 suggest a link between the presence of the integral 
operator (d/dx)^^ and the possibility of breaking the symmetry between the n 
and n—1 nodes of the Dn Dynkin diagram by choosing 7„ ^ 7n-i in the BAE. 
Therefore, in writing (16. ID we have omitted the integral operator {d/dx)~^ of 
(jS.ip because, in contrast to the Dn Dynkin diagrams, the Bn diagrams have no 
Z2 symmetry. 

The relevant solution iIj{x, E, g) to (|6.ip has the asymptotic a; ^ 00 behaviour 

^(x,ii;,g)~AAx-^ *^/2exp(^-^— , {M>K/{h''-K)) (6.2) 

on the positive real axis. The solutions (j3.22p 

Xi{x,E,g) , i = 0,l,...,2n-l (6.3) 
are instead characterised by the x — > behaviour 
y (x E ^) r^x^^+0(x^^+^^) {^i = 9i, fori<n-l. 

In (j6.4p the A's represent the 2n solutions of the indicial equation in table [2] with 
the ordering 

gi < gj < /i^/2 , Xi<\j , y i<j . (6.5) 

6.1 The ip-system and the Bn Bethe ansatz equations 

The Bn V'-system is 

= ^ a = l,...,n-l, (6.6) 

2 2 

^(n-l)^(n-l) ^^j^(n)^^(n)] (5,7) 
4 4 4 4 

Using the identity (I4.18P we can express V'^"'' with a>l in terms of ^/^'•^^ = ^. The 
result is 

= , a = l,2,...,n. (6.8) 
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Following the derivation in ^4.3^ define Ua = Yl'i=Q Aj, 

Qt\E,^) = Q[oi.,a-l](f^'^'g) ' Qt\E,^) = Q[oi.,a-2,a](^'^'g) (6-9) 

where and a-2a] defined as in (j4.29p and 

^t\^,E,^) = u;'"'^Qi^\E,g)W'^[xo,...,Xa-i] + 

^Ka^+i~X.-i)Qi-)(^E, g)W-[xo, Xa-2, Xa] + ■ ■ ■ ■ 

(6.10) 

Using (jG.lOp in (j6.6p and identifying the leading contributions about x=0 gives 
Q('^+i)(^)Q('^-i)(E) =a;^(^'^-^«-i)Q5°i(^)Qi"^(^)-'^^^^""'~^"^<3i"H^)Oil(^) 

2 2 2 2 

(6.11) 

with Q(°)(S) = 1 and 

2 Vl J 2 _ _Q-27a+7a_i+7a + i /g 

2 2 

In (j6.12p . a = 1, . . . , re — 1 and 



7a = a 



Aj + aw , = 1, 2,..., re, (6.13) 



where a = and u is still to be fixed. Plugging (jG.lOp into (j6.7p leads to 
Q[_Y'^(E)Q(r'^(^) = u;3(^"-^'-i)Q^1(^)Q(")(£;)-a;3(^"-i-^")g(")(^)Q^1 (S) 

4 4 4 4 4 4 

(6.14) 

and 

2 2 

The boundary condition jn+i='yn-i fixes u = — /i^/2, the {^a} {7a} relation 
in tableland allows (j6.12p and (j6.15p to be recast in the form (|2.ip . 

We have checked the consistency of the re=2 case both numerically and an- 
alytically. The relation with the sine-Gordon model briefly mentioned at the 
beginning of ^and the analysis of §6.21 and §7. 1.21 lend extra analytic support to 
the proposal. 
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6.2 Example 3: Bi 

This is a singular limit for the analytic BAE study in [37]. It however suggests 
the equivalence of the K=l case of Bi to the K=2 case of Ai in the integrable 
system. Indeed, the differential equation is second order in the former case and 
it can be written in the more-standard form 



Performing a Liouville transformation 



(6.16) 

il;{x,E,go)^i;{x,E,go)exp(^-^j (6.17) 

we find 

Equation (j6.18p coincides with the equations studied by Sergei Lukyanov [38] 
which are related to the Ai lattice models with integer spin 

j = K , K = 1,2,3,... . (6.19) 

The cases go=0 with M=K can be solved in closed form. After a shift x 
X + E, the simplest case K=l becomes 

which has general solution 

y{x) = ci e-'^H,! (^) + C2 e-^y^ I_i (^) (6.21) 

where and I are respectively the Hermite and the Bessel functions. Since 

F^(x) ~ 2'^x^(l + 0(l/x)) (6.22) 
as Re(x) > 0, oo, while 

/,(x)~-^(l + 0(l/x)), (6.23) 
V 27rx 

the most subdominant solution at large x on the positive real axis is 

^{x) = e-4H_i (^) . (6.24) 
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Figure [Tj shows the first five complex-conjugate pairs of zeros of ^(— e^). This 
2-string pattern is typical of ^i-related spin-1 integrable models. The exact 
eigenvalues are reported in table [5] above. There is also good agreement between 
the position of the first pairs of zeros shown in figured] and the WKB asymptotic 
prediction of appendix [Al 
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Figure 1: 2-strings in the Bi model at gQ=0, M=K=1. Contour plot of 1/(1 + 
\H_i( — 7^e^)|) in the complex 0-plane. 

2 V 2 

7 The Cn models 

The analytic and numerical results of ^ support the conjectured link between 
the Dn BAE ([23]) and equation (HH). At g={0, 1, . . . , n - 1}, the ODE is 



•2n-l 



d 



:j - VPk{x,E) ( — ) ^/PKix^E) ] ^(x, ^) = . 



(7.1) 



In this section, we start from (j7.ip and consider a second duality relation 

{D-n)K X {D-n)L {Cn)-K/2 X {Cn)~L/2 



{D-n)K+L 



{Cn)-K/2-L/2 



(7.2) 



discussed by Hornfeck [41j. Motivated by the results of ^4.11 on the analogy 
between the An <-> A-n spectral duality and the An-icelated duality in ^4ij, we 
change n — > —n and K —2K in (|7.1|) : 



-1 d 



-2n-l 



d 



^^-2n- + VP-2K{x,E) (^^^ ^P.2K{x,E)j^p{x,E)=0, (7.3) 

(7.4) 



where 



M(2n+2) \ —2K 

P-2k{x,E) = [x^^ -E' 
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Replacing 

ij{x,E) ^ (^{P^2K{x,E))-h (^A^ {P_2K{x,E)y^'^ Hx,E), (7.5) 

multiplying by ^^'2n+i and noticing that 

{P^2k{x,E))-'2^Pk{x,E) = (x'^'^^/^-i?)^, (7.6) 

where =n + 1 is the dual Coxeter number of C^, we can write the resulting 
equation as 



-1 



\dx^n+^ ^^^"'^^UJ PK{x.E)y{x,E)=Q. (7.7) 

Equation (|7.7p is our C„ candidate at g=go={0, 1,2,...}. Further, adapting the 
discussion of ^that led to the full Bn equation, and noting the similarity between 
the pseudo-differential operators (j3.18p . ()3.19p and ()3.20p . we replace 

ffn+l . / d \ / d \ 

and the final C„ proposal becomes (|3.2ip : 

[Dnig^] (-^) Dn{g) - Pk{x, E) (j^ Pk{x, E)] V(x, E,g)=0. (7.9) 



The relevant solution of ()7.9p has the asymptotic x ^ oo behaviour 

^{x,E,g)^Mx~^^exp(^-j^j , (MyK/ih-'-K)) (7.10) 

on the positive real axis. The solutions (|3.22p 

Xi{x,E,g) , i = 0,l,...,2n + l (7.11) 
are characterised by the x ^ behaviour 

{Xi= gi , for i < n - 1, 

A„ = n, 
Xi = 2n- g2n-i , for i > n . 

(7.12) 

In (j7.12p the A's represent the 2n + l roots of the indicial equation in table [2] with 
the ordering 

Qi < gj <n , Xi < Xj , M i < i . (7.13) 
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7.1 The ^/'-system and the C„ Bethe ansatz equations 

We now deduce the Cn Bethe ansatz equations from the proposed C„ ■i/'-system: 
VF[V'i1,V'i''^] = a = l,2,...,n-2 , 

4 4 

l^fViV'^V'l""'^] = , (7.14) 

4 4 4 4 

2 2 

Using the Jacobi identity we find 

0W =VF(")[V'i^,V'3^---,V'ii^] , a= l,2,...,n-l , (7.15) 

4 4 4 

^^^''iVi''^ = W^(")[^ 1^,-03^,..., V'li^] , (7.16) 

4 4 4 4 4 

where the functions (f)^'^\x , E , g) satisfy the fohowing system of functional rela- 
tions 

W[(l)'i\ , (A^"^] = i^^^'^+i) , a = 1, . . . , n - 1 , (7.17) 

4 4 

and 

W[<P^^l , (/.^r^] = 0("-i)(V'("))2 . (7.18) 

4 4 

Prom (I7.18P we see that 



(x, i?, g) = . (7.19) 

Now set 

(x,i?,g) = a;'=-"Qi'^)(i?,g)H^(«)[xo,...,Xa-i] + 

^fcK+,-A..,)QW (ij;, g)vFW [xo, . . . , Xa-2, Xa] + . . . (7.20) 

with a = l,2,...,n. The orders of the first and the second term in (j7.20p are given 

respectively by aa — a{a — l)/2 and aa+i — ^a-i —a{a — l)/2, where = Yli=o ■^i- 
Assuming 

^i") {x, E, g) = c^'^^x^qI") {E, g) + o.'^'^x'^qI") {E,g) + ... (7.21) 
where P < a, we can make the following identifications 

/3 = a„/2 - n(n - l)/4 , a = /3 + (A, - A„_i) (7.22) 

and 

Q(0)(S) = Q(0)(E) = 1, (7.23) 

Q(")(S) = Q^^Ke) , a = l,2,...,n-l , (7.24) 

Q(")(i?) oc Q^^l{E) Q^r\E) , (7.25) 

4 4 

QW(i?) oc u;3(^"-^"-i)Q5_1(ii;) g(r^(i?) +u;i(^"-^-^")gl")(ii;) Q^"](ii;). 

4 4 4 4 

(7.26) 
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Using ()7.20p in (j7.17p and selecting the leading terms we find 

4 4 4 4 

(7.27) 

With the identifications ([7^ and fr25]) this leads to 

rr ^^"^^^^ % = -j7f (^-^«-) , (7.28) 

where a = 1, 2, . . . , n - 1 and a = . Using ([7:20]) with a=n in (HH]), ({7:26]) 
and ()7.25p gives instead 



Q("-i)(^) =c^|(^"-^"-i)Q("|(^)Q(^")(^)-cji(^"-i-^")Qf)(^)Q(_"|(E); (7.29) 

2 2 2 2 

which leads to 



(n-l)( (n), 

, = . (7.30) 

2 

Finally, with the identification (|7.12p and the choice in table[3]for the {ga} ^ {la} 
relation, equations ()7.28p and (I7.30|) can be assembled into the universal form 
(O). 



7.1.1 Example 4: Ci 

The n=l case is again a singular limit of the analytic BAE, but it also suggests 
the similarity of this case to the Ai models [37j . The pseudo-differential equation 
is, however, not second order but instead third order 



hr^-— :r + — -^^(^'^)hr PK{x,E)\^i:{x,E,gQ) = {) (7.31) 

\ ax-^ x'^ ax x'^ \dx J I 

where C = go{go — 2). It is nevertheless easy to check that (|7.3ip is solved by the 
product of two functions satisfying second order ODEs: 

i^{x,E,go) = x-{x,E,go)x+{x,E,go) (7.32) 
where x± originates from a single function x as follows 

X±{x,E,go) = x{uJ^'^^x,n^'/^E,go). (7.33) 
Since we assume that x satisfies the standard ODE associated with Ai , 

^-^PK{x,E)-^^x{x,E,go) = 0, (7.34) 
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the functions x± satisfy the following: 



^^2^^^^'^'9o) = y±^PK{x,E) + ^ j x±{x,E,go) . (7.35) 
Starting from (j7.33p and differentiating ijj = iJj{x,E,gQ) three times we find 

dx dx dx ' dx"^ dx dx 2x'^ ' 

d^i^ . D . dx- dx+ . C dij) C 

-r^ = iPk [X+-; + —-; 3 V' ■ (7.36) 

dx'^ dx dx dx x^ 

We notice that 

ax V ax dx 



and therefore 



dx dx \dx 



Inserting (|7.38p into (j7.36p we finally arrive at equation (j7.3ip . 
7.1.2 Example 5: C2 ~ B2 

The ODEs for B2 and for C2 are both deduced from the D„-type ODE, but 
through completely different routes. It is thus a good test to check the equivalence 
of these two cases. We start with the i32-related ODE at g={0, 1} 

dV ^ p di' dPK , „ 

The ODE associated with the second node of B2, which is nothing but the first 
node of C2, would be satisfied by the following function 

t/,(2) = iy[V'_i,V'i] = [0,1] , (7.40) 

2 2 

where 

We then easily evaluate the derivatives of ^/j^^^: 

#(2) ^2^(2) ^3^(2) 

^ = [0,2], ^ = [1,2] + [0,3], ^ = 2[1,3]+P;,^(2) 

= 2[2,3]+P^^, ^ = Pk^-2P^[1,2]. (7.42) 
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Finally, noticing that 



211,2] = 2 (^-j [l,31 = ^-(^-j Prf^' (7.43) 

we have obtained a closed form ODE for ^(^^i 

' PKix,E)^ip(^\x,E)=0. (7.44) 

This equation is exactly the g={0, 1} equation associated to the first node of C2. 
This verifies the consistency of our proposal. 

More generally, the correspondence between the B2 parameters E={goy 9i} 
and the C2 parameters E={9o,9i} is established by 

2^0 = 50 + ffi - 1 , 

2gi = 50-91 + 3. (7.45) 

Finally we have directly checked the C2=B2 V'-system at g={0, 1} by verifying 
that both sides of 

W[i;^y, , Vl^^] = i^^^l V'l^^ , (7.46) 
44 44 

satisfy the same 15th order ODE. 

Unfortunately we do not currently have NLIEs for the Bn/Cn models and the 
numerical checks presented below are instead based on an approximate solution 
similar to that used by Voros in [5]. For the B2 BAEs for the case g={0, 1}, K=\ 
and Af=2/3, we started from a perfect-string estimate for the first 1000 roots as 
follows: 

^. = (4y|^j)^ node (7.47) 

Ej = e±*f (4V^^(j - h)-^ 2nd node. (7.48) 

1(3) " 

We then solved the BAEs recursively using the Newton-Raphson method on 
the first 20 roots, keeping the remaining roots fixed. Table [7] compares the lowest 
roots thus obtained with the results from the solution of the (pseudo-)differential 
equations. The relatively low accuracy in comparison with tables [H [5] and [6] is 
most likely to be a consequence of the slow convergence rate of the algorithm 
used to solve the Bethe ansatz equations, and in particular the systematic errors 
introduced by fixing the higher levels {Ej, j > 20) to their perfect-string values. 
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BA numerics 




ODE numerics 




1st node 


2nd node 


1st node 


2nd node 


6.28405 


6.8368 ± 5.8640i 


6.28390 


6.8365 ± 5.8637i 


13.2379 


18.216 ± 14.265i 


13.2376 


18.214 ± 14.264i 


21.6307 


30.996 ± 23.645i 


21.6303 


30.992 ± 23.642i 


30.5039 


44.747 ± 33.707i 


30.5034 


44.739 ± 33.700i 


39.8617 


59.254 ± 44.304i 


39.8613 


59.240 ± 44.292i 



Table 7: C2=B2'- comparison of Bethe ansatz results with the numerical solution of the B2 and 
C2 equations using the algorithm described in appendix|B] {B2 node convention with g={0, 1}, 
K=l and M=2/3.) 



8 Conclusions 

There are many aspects about the correspondence between integrable models and 
the spectral theory of ordinary differential equations that we would like to explore 
and understand at a deeper level. Given the applications of the Bethe ansatz to 
the study of QCD in its leading-logarithm approximation [50] and to the study 
of anomalous dimensions of composite operators in Yang- Mills theories [511 152j . 
the extension of the correspondence to lattice models is certainly desirable from 
a physical point of view. On the other hand the mathematical structures arising 
from the generalisation of the correspondence to other conformal field theories 
with extended symmetries has the potential to link areas of modern and classical 
mathematics in an elegant way. In this paper our results were obtained very 
much on a case-by-case basis, but we already saw the emergence of interesting 
mathematical objects: the V'-systems, negative-dimension dualities, and the for- 
mal similarity with the Miura-opers studied both in the classical work jl9] (see 
also [53l [Ml ES]) and more recently in [HI [H HSl [HI HE] . 

It would be very interesting to generalise equations (|3.181[3T2T|) to encompass 
the excited states of the integrable models; to date this has been completed for the 
K=l case of Ai \56\ I57j. More challenging, but also extremely interesting, would 
be to extend the correspondence to perturbed conformal field theories defined on 
a cylinder, both for the ground state [58] and for excited states [59| [60]. 

Finally, even remaining inside the current setup, the ODE/IM correspondence 
has already had an impact on condensed matter physics: it has been applied to 
interacting Bose liquids [6IJ, the single electron box [62j and quantum dots [63] . 
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A Root strings from the complex WKB method 

As mentioned in ^ a characteristic feature of the 'fused' models with K>1 is 
that the asymptotic roots are not necessarily real, but rather are grouped into 
so-called 'strings' of complex roots. Here we show how this behaviour can be 
recovered from a treatment of the ODE using the complex WKB method. We 
shall restrict the analysis to the Ai case and for simplicity set g={0, 1}. The 
ODE is the n=2 case of KT\i : 

[-^ + PK{z)^i^{z,E)=0, PK{z) = {z'^'/''-Ef (A.l) 

and the boundary condition determining the eigenvalues is that there should be 
a solution which decays as z ^ oo on the positive real axis, and is simultaneously 
zero at z=0. (In this section we use z instead of x to emphasise that it must be 
considered in the complex plane.) 

Before giving the WKB treatment, we summarise the expectations from inte- 
grable models. A useful discussion of the K=2 case, for generic twist parameter, 
is in appendix 2 of [6l]. For a finite lattice of size A^, the integral equation 
derived in that paper leads to the following asymptotic condition on the Bethe 
ansatz roots 9^ 

N In f tanh f \ = ±Nm =F {2j- 1 H —)TTi - In 2 , (A.2) 

\ K'^'^ J J vr-27 2 

where 1 < j < N/2 + 1. This equation describes the approximate position of 
complex-conjugate pairs of roots {^Ji^j*^}! with 9j and 6^ lying in the upper 
and lower halves of the complex plane. The relationship between 7, (j) in (IA.2P 
and M and g={5'o, 1 - Qo] is 

M+l TT 1 

^ = ^^- = 2^^' ^ = 2"^°- (^-'^ 

In the continuum limit, N ^ 00 and the numbers e ^ with 1 < j ^ A^/2 tend to 
zero. Shifting 9^ appropriately, equation (IA.2I) becomes 

e^'^'f = ±{2j - i - go)7ri - ^ In 2 , j = 1, 2, . . . . (A.4) 
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The second term causes deviations from the string pattern, even within this 
asymptotic approximation. Set 9^ = ±z7r/(4//) + ; then 

e2M/5^ = (2j-i-5o)vr±^ln2 (A.5) 
or, taking logs and expanding for large j, 

2fi(3f = ln(2j -I- 5o)vr ± i + . . . . (A.6) 

^ 2 2(2j - i - go)ir 

Equation ()A.6p exhibits the asymptotic deviations from the perfect string config- 
urations, which are only recovered in the limit j — > oo. This qualitative pattern 
matches the results of table [5] and figure [H and is illustrated in figure [2j To 
compare with the Ai , K=2 results in ^4.31 one has to identify 

^^2^.9+^^2^.e-^ = {^,-\E2,-iY , ^l-\E2,Y] (A.7) 



m 



foil . 



^ O O O— -0---0--0--G-G^O 

1 .5 - 



Im 1: 
0.5- 



0; 1 1.5 2 2.5 3 3.5 4 4.5 
Re 



-0.5- 



-1.5-__ 



"D o o — o— o-o-e-o-©-eH3-oesoQs 



Figure 2: Two-strings from the asymptotics of the NLIE |64j . 

These features (and their generalisations to higher K) can all be recovered 
from a complex WKB treatment of (|A.ip . The practicalities of this method have 
been very clearly explained by Heading [651 [66], s-nd we shall, more-or-less, follow 
his notations, which we now summarise. The behaviour of (jA.ip is controlled 
by Pk{z) = (z^M/Zs: _ ^ ^j^g possibly-complex zeros of this function are 
called turning (or transition) points. For each turning point zq, a branch cut 
is inserted joining zq to infinity; and for 2M/K ^ Z, another is added starting 
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from the origin along the negative-imaginary axis. After this, taking arg {Pk{z)) 
to tend to zero along the positive real z axis renders {Pk{z))^^'^ and (Pft'(z))^/^ 
uniquely defined. Prom each turning point zq there is a set of lines 



[ dt^Pxit) = (A.^ 

J Zt) 



called Stokes lines, and another set 

5Re / dt^/PK{t) = (A.9) 

J Zo 

called anti-Stokes lineS. A K^^ order turning point has K+2 Stokes lines, and 
K+2 anti-Stokes lines, emanating from it. Leading-order WKB solutions are 
written as 

{zQ,z) = (Pi^(z))-i/^exp (^jyt^PKit)^ . (A.IO) 

If 3?e j^^ di-sj PK{t) > 0, a subscript d is added, signifying that {zq, z)d is a domi- 
nant solution; (z, zq) is then subdominant, and is written (z, zo)s- On anti-Stokes 
lines, dominant and subdominant solutions swap role^. 

Formal WKB solutions a{z,zo) + (3{zo,z) are only valid in restricted do- 
mains: the Stokes phenomenon means that the coefficient of the subdominant 
term changes by an amount proportional to the dominant term when a Stokes line 
is crossed, so that a {z, zo)ii + P {zq, z)s is replaced by a {z, zo)ci + {P + Ta) {zq, z)s , 
where the Stokes multiplier T is a constant characteristic of the given Stokes line 
and crossing direction. (Note, since the discontinuity always occurs in the coef- 
ficient of a subdominant term, there is no contradiction with the fact that the 
WKB solution provides an asymptotic approximation to the exact - continuous - 
solutions of the original equation.) For a Stokes line emanating from a K^^ order 
turning point, traversed in a positive (anticlockwise) sense, the Stokes multiplier 
is [67] 

T = 2i cos{tt/{K+2)) . (A.ll) 

This can alternatively be recovered via the £'=0, /=0, M=K/2 value of the 
function C{E,l) discussed in [8], equation (3.10). For a clockwise traverse, T is 
replaced by — T. 

The key feature of (|A.ip . leading to the formation of strings, is that the 
turning points which control the WKB eigenvalues all have order K. If a solution 



^Here we follow the conventions of Heading and, for example, Berry, but beware that other's 
conventions are exactly the reverse. Note also that our Pk{z) is minus the function q{z) used 
by Heading. 

" Note, in contrast to elsewhere in this paper, the terms dominant and subdominant are used 
here with respect to the WKB expansion parameter, which for brevity we have set equal to 1. 
The parameter, a factor of in front of the derivative term in ()A.1|I . can easily be restored if 
desired. In Heading's notations, e corresponds to k~^. 
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which decays along the positive real axis is traced in to a turning point, there 
are K different directions away from that turning point along which the two 
WKB solutions (zq, z) and {z, zq) might come into balance and have a chance to 
cancel; if the continuation of any of these directions can be arranged to hit the 
origin, then a (WKB) eigenvalue is possible. When the Stokes multipliers have 
unit modulus, these directions are precisely the anti-Stokes lines; departures from 
unit modulus cause small shifts, and lead to deviations from the string pattern. 



Im 




-5- 



5 10 
Re 



15 



Figure 3: Initial directions of Stokes lines (continuous) and anti-Stokes lines (dot- 
ted) from a 4*^ order turning point. 

An example should make this clear. Figure [3] shows the set of Stokes and anti- 
Stokes directions away from a turning point of (1A.1[1 with M=3, X=4. Between 
each pair of anti-Stokes lines there is a Stokes sector, which is further divided into 
two segments by the Stokes line which lies exactly in the middle of the Stokes 
sector. The first few of these segments are labelled 1 . . . 5 in the figure; note that 
segments 1 and 2 make up one Stokes sector, and segments 3 and 4 another. For 
the figure, arg(ii^) was chosen so that the origin lies in the union of segments 4 and 
5, and the positive real axis in the union of the continuations of segments 1 and 2. 
Across the (dotted) anti-Stokes lines, dominant and subdominant solutions swap 
over, while across the (continuous) Stokes lines, the coefficients of subdominant 
terms may change by the Stokes phenomenon. Figure [4] shows the continuation 
of the Stokes directions shown in figure [3] into the full complex plane, justifying 
the claim that the real axis lies in the continuation of segments 1 and 2. 

Our aim is to find an approximate solution which decays as z — > oo on the 
positive real axis, and is also zero at the origin. Since the positive real axis lies 
in the continuations of segments 1 and 2, we start with a solution subdominant 
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Figure 4: The full Stokes and anti-Stokes lines for the situation shown in figure [3j 

in segments 1 and 2, and continue it clockwise around the turning point to find 
the behaviour in the region of the origin. There is no Stokes phenomenon passing 
between 1 to 2 since the dominant term is absent in both of these segments, but 
from then on the phenomenon must be taken into account each time a Stokes line 
is crossed. With zq the location of the turning point, and {z, zq) the approximate 
sub dominant solution in segments 1 and 2, we have: 

1,2: {z,zo)s 

3 : iz,zo)d 

4 : {z,zo)d-'2icos{TT/6){zo,z)s 

5 : {z,zo)s-2icos{TT/G){zo,z)d 

Near to the anti-Stokes line between segments 4 and 5, the WKB solutions 
{z, Zq) and 2i cos(7r/6) (zq, z) are of similar magnitude, and the (WKB) condition 
yielding the eigenvalues is found by demanding that they cancel exactly at z=^: 



(zo,0) - 2icos(7r/6)(0,zo) = 0. 



(A.12) 



That is 




(A.13) 




(A.14) 



Remembering that zq=E^/^'^^\ the integral can be evaluated to 





{-1)^/'^k{2M/K, 2/K) Ef" 



(A.15) 
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where ^ = K{M+l)/{2M) and 



r(i + l/a)r(l + 1/b) 

r(i + i/a + i/b) ■ (^-''^ 

Equations (|A.14I — [A.16|) give the WKB prediction at K=A. Had we instead per- 
formed the calculation for K=2, the logarithm on the RHS of ()A.14p would have 
been replaced by — ln(2 cos(7r/4)) = — ^ln2, and ()A.4j) would have been repro- 
duced. More generally, the Stokes multipliers solve a stationary T-system, and 
the constant determining the deviation of the j^^ root in a iiT-string is found to 
be 

ln(sin(^i/(K + 2))/sin(7r(i + l)/(ir + 2))), j = 1, 2, . . . , i^. (A.17) 

The reason why only K out of the K+2 anti-Stokes directions allow for zeros of 
the wavefunction is that on the two anti-Stokes lines next to the sector where the 
initial subdominant solution is defined (the anti-Stokes lines bordering the union 
of segments 1 and 2 in the example) there is only a single 'pure' WKB solution, 
and thus there is no chance to get the cancellation between two WKB solutions 
which leads to the wavefunction zeros near the other anti-Stokes directions. 

It would be interesting to reproduce these results from a study of the relevant 
nonlinear integral equations in the integrable models, but we will leave a more 
detailed study of such issues for future work. 

B The numerical algorithm and the dual formulation 
of the boundary problem 

In this appendix we describe the numerical methods used to test our conjec- 
tures. We solved the pseudo-differential equations using an iterative power-series 
method, which we describe below. On the integrable model side, there is a well- 
known method to rewrite an infinite set of Bethe ansatz equations such as (12. ip 
into a finite set of nonlinear integral equations (NLIEs) |68j . We applied this to 
the K=l simply-laced cases A„_i and Z?„. 

B.l The x-functions and the generahsed Cheng algorithm. 

The special solution tp{x, E, g) and its derivatives can be written in terms of the 
X-functions ()3.22p and the spectral determinants as 



^^{x, E, g) = Q W {E, g)^Xo(x, E, g) + Q W {E, g)^Xi {x,E,g) + ... . 

(B.l) 
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By solving for Q[g]^(^,g): 

In order to calculate the spectral determinant from ()B.2p we must first find the 
X-solutions by numerically solving the pseudo-differential equations. This can be 
done, with very high precision, using a generalisation of the iteration method 
introduced by Cheng many years ago [20] and more recently used in [69]. In the 
following we shall define the iterative solution for each Lie algebra on a case-by- 
case basis. 

An-i'- We begin by defining a solution for (j3.18p 

Dn{g)x{x, E, g) = {-iTPKix, E)x{x, E, g) , (B.3) 
and its adjoint equation 

Dnig^)xHx,E,g^) = PK{x,E)xHx,E,g^) , (B.4) 

where g={gn- i, ■ ■ -,91,90} and g^={n-l-go,n-l-gi, . . . ,n-l-gn^i}. 

Equation ()B.4p was defined, as usual, by applying the rule {dP/dx^)'^ = 
{-IfldP/dxP). Notice also that (iRi]) is not the -related ODE, which 

is instead obtained by simply replacing g with g^^ in (|B.3p . 

Consider equation (|B.3p . The first step is to define a linear operator 



X 



p+n 



L^(xP) = . (B.5) 

Wb=Q ip + n- 9b) 

such that for any polynomial V{x) of x 

D^{g)L^{V{x))=V{x) . (B.6) 
Now it is easy to check that 

Xa{x, E, g) = x'^- + i-lTL^ {Pk{x, E)xa{x, E, g)) (B.7) 
defines a solution to (|B.3p solvable by iteration. Setting gt = n — 1 — gn-i~a, 

x\{x, E, gt) = + L^t {Pk{x, E)xi{x, E, gt)) , (B.8) 
instead defines an iterative solution to the adjoint equation ()B.4p . 

Dn- The appropriate linear operator for the self-adjoint equation (j3.19p 

Dn{g^) (^) Dn{g)x{x, E, g) = y'PKix,E) (^^) ^Pk{x,E)x{x, E, g) 

(B.9) 
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IS 

and the corresponding iterative solution is 

Xa{x, E, g) = x^" + [y^PKix, E) ^/Pk{x, E) Xa{x, E, g)) . (B.ll) 
The relation between the A's and the g's is given in (|5.4[) . 

Bn- The iterative solution of the Bn equation ()3.20p 

D.n{g^)Dn{g)x{x,E,g) = -y/PK{x,E) (^^^ Pk{x , E)x{x , E , g) (B.12) 
and its adjoint 

D4g^)Dn{g)xHx,E,g) = VPk{x,E) (^^^ /FV(^X^(x,i?,g) (B.13) 
are, respectively, 

Xa(a;, ^, g) = x^- - (/fM^) VPk{x,E) xa{x, E, g)) (B.14) 
1 

xl(x, ^, g) = x^- + (^y/PKix,E) (^^^ VPk{x,E) xi{x, E, g)) , (B.15) 



and 



with 



The relation between the A's and the g's is given in (16. 4p . 
C„: The x-solutions to the self-adjoint equation (j3.2ip 



^n(g^) (^^ J ^n(g) Xa(x, ii^, g) = Pk{x, E) — J Pk(x, E) Xaix, E, g) 

(B.17) 

satisfy 



dx 



Xa{x,E,g)=x^^ +L'^(PKix,E){^] PKix,E)xa{x,E,g)] (B.18) 



with 



lC'(xP) = , (B.19) 

n^:o'(p+i+A.) 



and the relation between the A's and the g^s given in (|7.12p . 
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B.2 Dual formulation of the boundary problem 

From ()B.2p , with the numerical estimates for the ^-functions and their derivatives 
at large values of x, replacing ip{x,E,g) by its asymptotic behaviour and by 
varying E one has in principle access to Q^{E,g) and in particular to its zeros. 
This process is certainly possible but it is tedious and the CPU time increases 
at least quadratically with the order of the equation. Surprisingly, there is a 
short-cut that makes the algorithm essentially order-independent. 
To see this, consider the An-i case 

CQ\;^{E,g) = W[^,Xi,X2,---,Xn-i] , (B.20) 

where 

n-l 

C = VF[X0,Xi,X2,...,Xn-i] = (-ir*["/'] n (^*-^^)- (B-21) 

i=Oj=i+l 

Expand ()B.20p with respect to the first column: 

CQ\SiE,g) = J2(-^r^P^' (B.22) 

p=0 



where, borrowing the short-hand notation of ^5.31 



Wp{E,g) = [0,l,...,p-l,p+l,...,n-l] 



det 



(B.23) 



with X = (xi; X2; • • • , Xn~i) ■ At large x on the positive real axis 



— V'(x) ~ (-If exp(-x*'^+V(M+l)) , {M > K/{h^-K)) . 

^ (B.24) 
Therefore, if for all functions Wp{E,g) we have 

^p(e,,g) = o(x(^-")^^/2g^p(^M+i/(^^^))^ ^ p = 0,l,...,n-l (B.25) 
for some E=ei, then from ()B.22p 

Q[o/(ei,g) = . (B.26) 

This shows that {cj} C {E^^^}. We shall now argue that the set {cj} is not empty 
and further prove that 

{ej = {^f^}. (B.27) 
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To do this consider Wn-i{E, g) = [0, 1, 3, . . . , n — 2] . By repeated differentiation 
and substitution of the ODE, one can prove as in [TT] that 



n-l 



«;„_i(ii;,g) = ((-irt[("-i)/2l II {g,-g,))xi_,{x,E,g^) (B.28) 

i=l,j=i+l 

i.e. Wn-i satisfies the adjoint ^„_i-equation (|B.4p and thus can be determined 
using (|B.8P with gl^_-^ = n - 1 - go. 
At large x , 

= SiE, g) + o(x(i-")*^/2)) exp (^) + . . . (B.29) 

where the dots in ()B.29P indicate terms that grow exponentiaUy at most as 

e-P i^os ) ^) . (B.30) 



Therefore, the condition ()B.25P for p=n—l is fulfilled for any value of E=ei such 
that S{ei,g) = 0. Since this a single condition imposed on a nontrivial function 
of E, we expect an infinite but countable set of solutions. 

For the next step, we want to argue that if the condition (|B.25P is satisfied 
for p=n—l it is simultaneously satisfied for p=0, 1, . . . ,n — 2. 

From (|B.29P we see that at 5(ej,g) = also 

^^x;„_i(e„g) =o(x^exp(x^^+V(M+l))) (B.31) 

for any finite complex number b. We now claim that all of the Wa can be written 
as a linear combination of Wn-i and its derivatives: 

Wp{x,E,g)= —^——nj^_,(^x,E,g), (B.32) 

a=0 

where the constants xj^'' depend only on g. It is very simple to see that (1B.32P 
is true in general. Write the adjoint ODE in the 'expanded' form 

d^ ^"-^ =^^d^ + ^^ + -- - + ^^(^' ^) J ^"-1 • (B-33) 
Differentiating once 

■^Wn-i = [0, 1, 2, . . . , n - 3, n - 1] = Wn-2 , (B.34) 
taking further derivatives and using (|B.33p the general result is 
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Solving ()B.35P recursively proves ()B.32p . 

Equation (|B.32p and equation (1B.31|1 together are equivalent to (|B.25p . To 
prove ()B.27p . substitute the asymptotic behaviours (|B.24p and 

^'n-l-. ~ £l^n-l{x) ~ 5(i?,g)x(l-"+2.)M/2 gg) 

into (|B.22p to see that 

CQ\^iE,g)=nS{E,g) . (B.37) 

The eigenvalues are obtained by solving xli^ii^, E, g^)=0 for fixed x, chosen large 
enough to capture the asymptotic behaviour of the solution while also ensuring the 
power series is reliable. This condition is different from the precise requirement 
(|B.3ip . but considering (|B.29P and ()B.30P it is clear that, provided x is very large, 
the error can be minimised and it selects the points on the complex i?-plane where 
the function S{E, g) is approximately zero. 

Similarly, for Bn the eigenvalues of the spectral determinant are found 

by solving (IB.ISP for X2n-i(^' s)- The situation is more complicated in the 
presence of the integral operators that appear in the Dn and the Cn models, 
and we have not yet fully completed the analysis. However, guided by the above 
and the results when the integral operator is absent, the obvious prescription to 
compute X2n-i{x-,E,g) from (IB.lip and X2n+i{x-,E,g) from (jB.lSh for Dn and 
Cn respectively works very well. 



B.3 NLIEs 

The Bethe ansatz equations of type An-i and Dn for K=l can be rephrased as 
a set of r nonlinear integral equations where r is the rank of the algebra p)8] (see 
also [TOlIMlEIlElElIIIlEaiB])- They are compactly written as 

/W(0) = -|^^-2^6oe^ (B.38) 



\JCi JC2 



6=1 

with hQ=mi sin ^ and 

sinh(7r^A;) 

27r ' sinh(f^(/iv^- l)A:)cosh(f^) 

The integration contours Ci and C2 run from —00 to 00 just below and just above 
the real axis respectively. The algebra dependence is encoded in the constant bo, 
the 'twists' {ja/oi}, and the deformed Cartan matrix 

r 2 a=b 
C^k) = . (B.40) 

I cosh(^) 
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Above, <ab> implies the nodes a and b of the Dynkin diagram are connected. 
The exphcit expressions for the inverse Cartan matrices can be found in |68j . 

By construction, the spectrum of the ordinary differential equation is encoded 
in the zeros 6j of 1+ exp(f^"'\9) that lie on the real axis of the complex-0 plane [T]. 
The eigenvalues of the relevant ODE are obtained via the relation 

Ej = e^^l^" . (B.41) 

The large-i? WKB-like behaviour can be obtained from the large-0 limit of 
these equations. In this limit the convolution terms in (|B.39P can be dropped, 
and imposing f^^\9j)=iTT{2j — 1) for integer j yields the semiclassical prediction 

2sin(^7r)m,(4"^r ~(2j+l--^7a)vr , j = 0,1,2,... (B.42) 
with a=2(/i- 
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